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Abstract 

We consider the problem of forecasting the regions at higher risk for 
newly introduced invasive species. Favourable and unfavourable regions 
may indeed not be known a priori, especially for exotic species whose hosts 
in native range and newly-colonised areas can be different. Assuming that 
the species is modelled by a logistic-like reaction-diffusion equation, we 
prove that the spatial arrangement of the favourable and unfavourable re- 
gions can theoretically be determined using only partial measurements of 
the population density: 1) a local "spatio-temporal" measurement, dur- 
ing a short time period and, 2) a "spatial" measurement in the whole 
region susceptible to colonisation. We then present a stochastic algorithm 
which is proved analytically, and then on several numerical examples, to 
be effective in deriving these regions. 

Keywords: reaction-diffusion ■ biological invasions • inverse problem • habi- 
tat configuration • Carleman estimates • simulated annealing 



1 Introduction 

Because of trade globalisation, a substantial increase in biological invasions has 
been observed over the last decades (e.g. Liebhold et al. PQ). These invasive 
species are, by definition |2J, likely to cause economic or environmental harm or 
harm to human health. Thus, it is a major concern to forecast, at the beginning 
of an invasion, the areas which will be more or less infested by the species. 

Because of their exotic nature, invading species generally face little compe- 
tition or predation. They are therefore well adapted to modelling via single- 
species models. 

Reaction-diffusion models have proved themselves to give good qualitative 
results regarding biological invasions (see the pioneering paper of Skellam [3], 
and the books @], [5] and [5] for review). 



The most widely used single-species reaction-diffusion model, in homoge- 
neous environments, is probably the Fisher-Kolmogorov [3 [8] model: 

u t = DAu + u(ji - 7it), t > 0, x G Q C R N , (1) 

where it = u(t, x) is the population density at time t and space position x, D is 
the diffusion coefficient, /i corresponds to the constant intrinsic growth rate, and 
^ is the environment's carrying capacity. Thus 7 measures the susceptibility to 
crowding effects. 

On the other hand, the environment is generally far from being homogeneous. 
The spreading speed of the invasion, as well as the hnal equilibrium attained by 
the population are in fact often highly dependent on these heterogeneities ([4], 
[TO] , [H]). A natural extension of |T]) to heterogeneous environments has 
been introduced by Shigesada, Kawasaki, Teramoto [12]: 

Ut = V(D(x)Vu) + u{(j,(x) - j(x)u), t > 0, leSlcK" (2) 

In this case, the diffusivity matrix D(x), and the coefficients n(x) and 7(2) 
depend on the space variable x, and can therefore include some effects of envi- 
ronmental heterogeneity. 

In this paper, we consider the simpler case where D{x) is assumed to be 
constant and isotropic and 7 is also assumed to be positive and constant: 

u t = DAu + u(u{x) - ju), i > 0, x e n C R N . (3) 

The regions where /1 is high correspond to favourable regions (high intrinsic 
growth rate and high environment carrying capacity) , whereas the regions with 
low values of /x are less favourable, or even unfavourable when \i < 0. In what 
follows, in order to obtain clearer biological interpretations of our results, we 
say that /i is a "habitat configuration" . 

With this type of model, many qualitative results have been established, 
especially regarding the influence of spatial heterogeneities of the environment 
on population persistence, and on the value of the equilibrium population den- 
sity ([4], [9], [13], [14], [15]). However, for a newly introduced species, like an 
invasive species at the beginning of its introduction, the regions where /1 is high 
or low may not be known a priori, particularly when the environment is very 
different from that of the species native range. 

In this paper, we propose a method of deriving the habitat configuration jj,, 
basing ourselves only on partial measurements of the population density at the 
beginning of the invasion process. In section [2] we begin by giving a precise 
mathematical formulation of our estimation problem. We then describe our 
main mathematical results, and we link them with ecological interpretations. 
These theoretical results form the basis of an algorithm that we propose, in 
section [3l for recovering the habitat configuration fi. In section [H we provide 
numerical examples illustrating our results. These results are further discussed 
in section El 
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2 Formulation of the problem and main results 



2.1 Model and hypotheses 

We assume that the population density u 7 is governed by the following parabolic 
equation: 

dfUj — DAu 7 + u 7 (p(x) — 7« 7 ), (>0,ieH, 
u J (t,x) = 0, t>0,x&dfl, (Pna) 
u~(0,x) = Ui{x) in fl, 

where fi is a bounded subdomain of M. d with boundary dfl. We will denote 
Q := (0, +oo) x fl and £ := (0, +oo) x dfl. 

The growth rate function p is a priori assumed to be bounded, and to take 
a known constant value outside a fixed compact subset S^i of fl: 

p £ M := {p £ L°°(fl), -M < p < M a.e., and p = m in fi\^i}, 

for some constants m, M , with M > 0; the notation "a.e." means "almost 
everywhere" , which is equivalent to "except on a set of zero measure" . 

The initial population density Ui{x) is assumed to be bounded (in C 2 (fl)), 
and bounded from below by a fixed positive constant in a fixed closed ball 
B e C fix, 01 small radius e: 

T> := {4> > 0, G C 2 (H), U\\ C 2 {n) < u-, <f> > Ui in B e , }, (4) 

for some positive constants ul and Ui . 

Absorbing (Dirichlet) boundary conditions are assumed. 

Remark 2.1 Absorbing boundary conditions mean that the individuals cross- 
ing the boundary immediately die. Such conditions can be ecologically relevant 
in numerous situations. For instance for many plant species, seacoasts are lethal 
and thus constitute this kind of boundaries. 

For technical reasons we have to introduce the subset fix, such that, in the 
interface between fix and ft, p. takes a known value m. This value is typi- 
cally negative, indicating that, near the lethal boundary, the environment is 
unfavourable. This assumption is not very restrictive since, in fact, fl\ can be 
chosen as close as we want to fl. 

For precise definitions of the functional spaces L 2 , L°° and C 2 as well as the 
other mathematical notations used throughout this paper, the reader can refer, 
e.g., to [TB]. 

2.2 Main question 

The main question that we presented at the end of the Introduction section 
can now be stated: for any time-span (to,t\), and any non-empty subset u> 
of fix, is it possible to estimate the function p,(x) in fl, basing ourselves only 
on measurements of u 7 (t, x) over (to,tx) x oj, and on a single measurement of 
u 7 (t,x) in the whole domain fl at a time T' = t ° J ^ tl ? 
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2.3 Estimating the habitat configuration 

Let p, be a function in M., and let v be the solution of the linear parabolic 
problem (P^o)- We define a functional G M , over K + x A4, by 

G M (7,M) = ll<9t" 7 -^Ili=((t ,t0x w ) ( 5 ) 
+ || An, (T', •) - Ati (T', •) || 2 L2(n) + K (T', •) - 5 (T', •) ||£ 2(f $i) 

where u 7 is the solution of (P Mi7 ). This functional G M quantifies the gap between 
u 7 and u on the set where m 7 has been measured. 

Theorem 2.2 TTie functions /i,jieM fcmj given, we have: 

C 

llM-Alli^no < =2^(0, /2), 
/or all jl £ M. and for some positive constant C = C(Cl,Cli,w,B £ , D,t ,ti, Ui/uj). 



The proof of this result is given in Appendix A. It bears on a Carleman-type 
estimate. 

Biological interpretation: This stability result means that, in the linear 
case corresponding to Malthusian populations (7 = 0), two different habitat 
configurations ji, ft cannot lead to close population densities Uo, v. Indeed, 
having population densities that are close to each other in the two situations, 
even on a very small region uj, during a small time period (to,ti), and in the 
whole space Cl at a single time T", would lead to small G M values, and therefore, 
from Theorem I2.2[ to close values of the growth rate coefficients /1 and jl. 

Theorem 12.21 implies the following uniqueness result: 

Corollary 2.3 If v is a solution of both (P M ,o) & n d (-P/i,o)> then fi — jl a.e. in 
ill, and therefore in CI. 

Biological interpretation: In the linear case (7 = 0), if two habitat con- 
figurations /j,, ji lead to identical population densities Uo, v, even on a very 
small region uj, during a small time period (to,ti), and in the whole space CI at 
a single time T", then these habitat configurations are identical. 

Next we have the following result: 

Theorem 2.4 We hav£\ |G M (0, jl) - G^, p)\ = C(W 3 ), as uj -> 0. 

The proof of this result is given in Appendix B. 

Biological interpretation: Assume that the habitat configuration n is not 
known, but that we have measurements of the population density it 7 , governed 
by the full nonlinear model ([3]) . Consider a configuration jl in A4 such that the 
population density v obtained as a solution of the linear model (-Po.ju) nas values 

'Two functions f(fi, fi, Ui,Ui,Ui, 7) and g(fi, jl, u% , Uj, W> 7); are written / = O(g) as g — > 
if there exists a constant K > 0, independent of fi, jl, Uj, ttj, ul and 7, such that |/| < K\g\ 
for g small enough. 
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close to those taken by the population density w 7 , in the sense that G M (7,/i) is 
close to 0. If the initial population density is far from the environment carrying 
capacity, then uj <C -, uj is small and, from Theorem 1 2. 4[ G M (0, fx) is also close 
to 0. Thus Theorem \2 . 21 implies that the habitat configuration jl is an accurate 
estimate of ji. In section [3j we propose an algorithm to obtain explicitly such 
estimates of /x. 

Remark 2.5 In fact, the term Oiuf) increases exponentially with time t\. 
Thus, obtaining accurate estimates of /i require, in practice, to work with small 
times i.e. at the beginning of the invasion. 

2.4 Forecasting the fate of the invading population 

The knowledge of an L 2 -estimate jl of [i enables us to give an estimate of the 
asymptotic behaviour of the solution it 7 of (P^ i7 ), as t — > +oo, and especially to 
know whether the population will become extinct or not. Indeed, as t — > +oo, it 
is known that (see e.g. [9], for a proof with another type of boundary condition) 
the solution Uj(t, x) of (P^ )7 ) converges to the unique nonnegative and bounded 
solution p-y of 

-DAp y = p y (fi(x) - 7p 7 ) in tt, (S^ y ) 

with pj — on <9f2. Moreover, p 7 = if and only if Xi[fj] > 0, where Ai[/z] is 
the smallest eigenvalue of the elliptic operator C : i— > —DAijj — fi(x)ip, with 
Dirichlet boundary conditions. On the other hand, if Xi\p] < 0, then p 7 = in 
Q (note that 7 does not appear in the definition of Ai). 
We have the following result. 

Proposition 2.6 Let us consider a sequence (p- n )n>o in M, such that fi n — » [i 
in L 2 (Q) as n — > 00. 

a) The solution p 1>n of the problem (5*^ nj7 ) converges to p 1 as n — > 00, 
uniformly in J7. 

b) Xi[fi n ] ~* Xi\p] as n — > +00. 

The proof of this result is classical and can be found in [HI [H] . 

Biological interpretation: Assume that the habitat configuration fj, is not 
known. We know that, in large times, the population density it 7 will tend to an 
unknown steady state p 7 (possibly 0, in case of extinction of the population). 
The part a) of the above proposition means that, if we know an accurate (L 2 -) 
estimate fi of //, then we can deduce an accurate estimate p 7 of the steady state 
p~ ll provided the coefficient 7 is known. Part b) shows that, even if 7 is not 
known, having an estimate jl of fi enables to obtain an estimate of Ai[/z], and 
therefore to forecast whether the species will survive or not. Indeed the sign 
of Xi[fi\ controls the fate of the invading species (persistence if Ai[/i] < and 
extinction if Ai[/z] > 0, see H31 [TJJ [T7] for more details) . 
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3 Simulated annealing algorithm 



Let (to,t±) be a fixed time interval, and u> C fii be fixed. We assume that we 
have measurements of the solution u~(t,x) of (P^, 7 ) over (to,t\) x u>, and of 
u<y( °t > x ) m ^- However, the function [i and the constant 7 are assumed to 
be unknown. Our objective is to build an algorithm for recovering fi. 

Remark 3.1 When the function u 7 is known, the computation of 6*^(7, •) does 
not require the knowledge of 7. 

The function /1 is assumed to belong to a known finite subset E of M. , equipped 
with a neighbourhood system. We build a sequence An of N elements of E with 
the following simulated annealing algorithm: 

n = 

Initialise Ao 
while n < N 

Choose randomly ci neighbour v of fi^^n 
if G^.v) < G^(7, An) 
An+1 *~ v 

else 

Choose randomly with an uniform law w G (0, 1) 

G M (T,An)-G M (-y,i^) 

if ui < e e <™> 

An+1 <- ^ 

else 

An+1 ^~ An 

endif 
endif 

n <— n + 1 
endwhile 

The sequence 0(n) (cooling schedule) is composed of real positive numbers, 
decreasing to 0. The simulated annealing algorithm gives a sequence An of 
elements of E. It is known (see e.g. [18]) that, for a cooling schedule Q(n) 
which converges sufficiently slowly to 0, this sequence converges in L 2 (fl) to a 
global minimiser fi of G M (7, .) in E (but see Remark 13.21) . 

Moreover, from Theorems 12.21 and 12.41 we have 

C C 

I|M-An||i2(n) < =2 G m(°' An) < =2 G u(7,An) +£l(«i), 

where £1 is a real- valued function such that £i(s) — > as s — > 0. Since fi <E E 
we obtain that, for n large enough, 

G u (7, An) < G M (7,/i), 



G 



and, from Appendix A 



—oGJj, fi) — * as Ui — » 0, 
the ratio Uijul being kept constant. We finally get: 

IIa* — Anlli 2 (o) ^ as 7*7 — > and n — > +oo, 

for a fixed ratio Ui/ul- Thus, for 7*7 small enough, and for n large enough, (l n is 
as close as we want to in the L 2 -sense. 

Remark 3.2 The cooling rate 0(n) leading to the exact optimal configuration 
with probability 1 decreases very slowly (logarithmically) and cannot be used 
in practice; see [TH] for a detailed discussion. Empirically, a good trade-off 
between quality of solutions and time required for computation is obtained with 
exponential cooling schedules of the type 8(n) = Oo x a n , with a < 1, first 
proposed by Kirkpatrick et al. [2Q]. Many other cooling schedules are possible, 
but too rapid cooling results in a system frozen into a state far from the optimal 
one. The starting temperature 9o should be chosen high enough to initially 
accept all changes p, n +i <— v, whatever the neighbour v. 

For this type of algorithm, there are no general rules for the choice of the 
stopping criterion (see [19]), which should be heuristically adapted to the con- 
sidered optimisation problem. 

4 Numerical computations 

In this section, in one-dimensional and two-dimensional cases, we check that the 
algorithm presented in section [3] can work in practice. 

In each of the four following examples, we fixed the sets fi, S7i and Ai and we 
defined a finite subset EcM equipped with a neighbourhood system. Then, 
for a fixed habitat configuration [i S E we computed, using a second-order finite 
elements method, the solution u(t, x) of (P M , 7 ), for D = 1, 7 = 0.1, t S (0,0.5), 
and for a given initial population density Uj. Then, we fixed to — 0.1, t\ = 0.4, 
and a subset ui C £li, and we stored the values of M( t "^" tl , x), for x G Q, and 
u(t, x) for (t, x) G (toj <i) xcj. Using only these values, we computed the sequence 
(/}„) of elements of E, defined by the simulated annealing algorithm of section[31 
the function fio being sampled arbitrarily, in a uniform law, among the elements 
ofE. 

In the following examples, we used 80 = 100, and 0(n) = 100 x 0.99" 
for the cooling schedule. Our stopping criterion was to have no change in the 
configuration /t„ during 500 iterations. 

The rigourous definitions of the sets E and of the associated neighbourhood 
systems which are used in the following examples can be found in Appendix C. 
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Figure 1: Example 1: a) exact habitat configuration fi; b) measurement of 
u 7 (0.25,a;) which was used for recovering (i; the domain u> is delimited by two 
black dots. The exact configuration /x was recovered after 1500 iterations of the 
algorithm. 



4.1 One-dimensional case 

Assume that n = (0,100), (li = [10,90], uj = [55,58], M = 2, m = -1 and 
Ui = 0.1(1 - x/100)sin(7rx/25) 2 . 

Example 1: the set E is composed of binary step functions taking only the 
values m and M. 

The function \i £ E and the measurement u 7 ( t "'^ tl , x) are depicted in Fig. 
[TJ (a) and (b). Two elements of E are said to be neighbours if they differ only 
on an interval of length 1 . 

The sequence (jj, n ) stabilised on the exact configuration /i after about N = 
1500 iterations. 

Example 2: the set E is composed of step functions which can take 21 
different values between m and M. 

The function [i 6 E and the measurement m 7 ( °"t 1 , x) are depicted in Fig. 
[2 (a) and (b). Two elements of E are said to be neighbours if they differ by 
(M — m)/20 on an interval of length 1. 

This time, the sequence (fin) stabilised on a configuration fi (Fig. (c)) after 
7500 iterations. The mean error in this case was J Q \^{ x ) ~~ fi(x)\dx = 0.05. 
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Figure 2: Example 2: a) exact habitat configuration /i; b) measurement of 
u 7 (0.25,a;) which was used for recovering /i; the domain w is delimited by two 
black dots; c) configuration fi := /U7500, obtained after 7500 iterations. 
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(a) (b) 

Figure 3: Example 3: a) exact habitat configuration /i. In the black regions, 
the depicted function takes the value 2; in the white regions, it takes the value 
— 1. b) measurement of w 7 (0.25, x) which was used for recovering fi; the domain 
u> is delimited by a black circle. The exact configuration \x was recovered after 
3000 iterations. 

4.2 Two-dimensional case 

Assume now that Q = (0, 20) x (0, 20), fli = [2, 18] x [2, 18], and that u> is the 
closed ball of centre (7; 7) and radius 3 . Assume also that M = 2 and m = — 1, 
and that the initial data is it, = O.lxy/400 sin(a;/4) 2 sin(y/4) 2 . 

Example 3: E is composed of binary functions which can only take the 
values m and M on each cell of a regular lattice. 

We fixed fi e E as in Fig. [3] (a). The measurement tt 7 ( * ^* 1 , x) is depicted 
in Fig. [3] (b) . Two elements of E are said to be neighbours if they differ only 
on one cell of the lattice. 

The sequence (/t n ) stabilised on the exact configuration \x after 3000 itera- 
tions. 

Example 4: E is composed of functions which can take 21 different values 
between m and M on each cell of a regular lattice. 

The configuration /j, and the measurement %( 2 ' x ) are depicted in Fig. 
0] (a) and (b). Two elements of E are said to be neighbours if they differ by 
(M — m)/20 on one cell of the lattice. 

The sequence (/}„) stabilised on a configuration fi after 14000 iterations (Fig. 
0J (c)). The mean error was -^j Jq° Jq° \h( x , V) — P>(x> y)\dxdy = 0.04. 
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Figure 4: Example 4: a) exact habitat configuration /u; b) measurement of 
u 7 (0.25,a;) which was used for recovering fj,; the domain ui is delimited by a 
black circle; c) approached configuration fiuooo- The darker the regions in (a), 
(b) and (c), the higher the values of the depicted functions. 
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5 Discussion and conclusion 



We have shown that, for an invasive species whose density is well modelled 
by a reaction-diffusion equation, the spatial arrangement of the favourable and 
unfavourable regions can be measured indirectly through the population density 
at the beginning of the invasion. More precisely, we considered a logistic-like 
reaction-diffusion model, and we placed ourselves under the assumption that the 
initial population density was far from the environment carrying capacity (it can 
be reasonably assumed at the beginning of an invasion). In such a situation, 
the position of the favourable and unfavourable regions, modelled through the 
intrinsic growth rate coefficient /z, may not be known a priori. This is especially 
true for exotic species whose hosts in native range and newly-colonised areas can 
be different. From our results, in the "ideal case" considered here, the position 
of these regions can be obtained through partial measurements of the population 
density. These partial measurements consist in two samples of the population 
density: 1) a "spatio-temporal" measurement, but very locally (in the small 
subset lj) and during a short time period and, 2) a "spatial" measurement in 
the whole region susceptible to colonisation (Q). 

The stochastic algorithm presented in section [3] shows explicitly how to re- 
construct the habitat arrangement /i from the above partial measurements of 
the population density. This algorithm was proved to be effective in both one- 
dimensional and two-dimensional cases, in section 01 through several numerical 
experiments. In examples 1 and 3, the algorithm converged to the exact habitat 
configuration. In examples 2 and 4, the sizes of the sets of possible habitat con- 
figurations were increased compared to examples 1 and 3. In those cases, the 
algorithm converged to configurations which were close to the exact ones. It is 
noteworthy that the spatial measurement in Q and the habitat arrangement /x 
can have very different shapes; therefore, /i cannot be straightforwardly deduced 
from this measurement. 

These results can be helpful in preventing biological invasions. Indeed, a sim- 
ple protocol, consisting of placing one trap in the invaded region, and recording 
the number of individuals captured by this trap over a short time-period (de- 
pending on the species characteristics), and performing a single survey of the 
number of individuals and their position in the whole considered region should 
allow, from our results, to detect the favourable areas, and to treat them pre- 
ventively. As we have emphasised in Proposition 12.61 a. the knowledge of an 
estimate of the habitat arrangement /i also allows us to forecast the final pop- 
ulation density, and therefore to detect the regions at higher risk, for instance 
in the case of harmful species. As recalled in Proposition 12. 61 b. having a good 
estimate of the habitat arrangement fi is also crucial to forecast the fate of the 
invasive species: persistence or extinction. 

On the other hand, we have to underline that our approach may not be 
adapted to some species, especially those which colonies are made of few in- 
dividuals. Indeed, the diffusion operator of our model can be obtained as the 
macroscopic limit of uncorrelated random walks. With such an operator, by the 
parabolic maximum principle [16j . it is known that, even with a compactly sup- 
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ported initial population density, the solution of our model is strictly positive 
everywhere on the domain as soon as t > 0. This means that the solution, and 
therefore the information, propagate with infinite speed, which is not realistic 
for discrete populations. This could induce a practical limitation of our method 
to a certain type of species, which are well modelled by continuous diffusion 
processes even at low densities (typically some insect or plant species, with high 
carrying capacity and growth rate). 

Note also that some of the mathematical tools used in this paper, and espe- 
cially Carleman estimates (see Appendix A), were initially not adapted to the 
nonlinear model considered here. Thus, we first considered, in Theorem 12.21 
the linear - or Malthusian - case. For populations whose density is far from the 
environment carrying capacity, the linear and the nonlinear problem have close 
solutions. In this situation, Theorem [2T4] extended the result of Thcorcm l2.2l to 
the nonlinear case of a logistic growth. 

The results of this paper could be immediately extended to the case of spa- 
tially varying functions "f(x). Another easy extension would be, for the "spatio- 
temporal" measurement, to use a partial boundary observation on a part T + 
of the domain boundary dft instead of sampling the population over a small 
domain u>. Using a new Carleman estimate (see we are indeed able to 
write a stability result for the coefficient /i, similar to that of Theorem 12.21 but 
with \\d n {d t u 7 ) - d n (d t v)\\l 2[ , t . )xr+) instead of ||<9 t u 7 - d t v\\ 2 L2([ta ti)XLj) in 
the definition of the functional G. 
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6 Appendices 

Let us introduce the following notations: for all t, t' G M, with t 1 > t, we denote 
Q\ = {t,t') x fl and E* = (t : t') x dfl. Throughout this section, with a slight 
abuse of notation, we designate by C any upper bounds in our computations, 
provided they only depend on the parameters f2, fii, to, B £ , D, to, ti, Ui/ul. 

6.1 Appendix A: proof of Theorem 12.21 
Carleman estimate 
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We recall here a Carleman-type estimate with a single observation. Let (3 
be a function in C 2 (51) such that 

1< (3 < 2 in fi, = 1 on dfi, min{\V/3(x)\,x G > and <9 n /3 < on dfl, 

where n denotes the outward unit normal to <9fi. For A > and t G (io,£i), we 
dehnc the following weight functions 

e \0(x) e 2A _ e A/3(x) 



= 71 7777- IT' x ) 



(t-*o)(*i-*) ' {t-t )(h-t) 

Let g be a solution of the parabolic problem 

<9 t <7 - L> Ag + a(a:)g = /(t, ir) in , 
g = 0on£ \ (P) 
g(0,x) = g (x) in Q, 

for some functions / G L 2 (Q^), and a, go G I/°°(f2). 
Then the following results are proved in [22] : 

Lemma 6.1 Lei q be a solution of (P). Then, there exist three positive constants 
Ao, Co and s > 1, depending only on CI, w, to awe? ti smc/i t/iat, for any A > Ao, 
£/ie next inequalities hold: 



a) ||M 1 (e-"g)||2 a(Qti) + ||M 2 (e- s "g)||i 2(Qti) + s\ 2 e - 2 *'V|Vg| 
+ S 3 A 4 / tl < Co [* 3 A 4 J* 1 e" 2s W + fl e~ 2 ^ (/ - aq) 



2 

2 



w/iere Mi and M 2 are defined by M\ip = —DAtp — s 2 \ 2 D\V (3\ 2 ip 2 ijj + s(d t r])^p, 
and M 2 ijj = d t ip + 2s\DipV/3.Vip + 2sX 2 D(p\S7/3\ 2 tp. Moreover, 



b) s~ 2 \~ 2 f Qll e-^ip-^dtq) 2 + (Ag) 2 ] + e- 2s ^\Vq\ 2 + s 2 \ 2 j Qt , e~ 2s ^\ 2 
s 2 A 2 fl J u e~ 2sr >^q 2 + S-U- 2 J Qtfl e- 2 ° s " (/ - aqf 



Stability estimate with one observation 

Let jU, jl G M. We consider the solutions v and v of the linear problems 
(Pfj,,o) an d (Pp,,o), respectively. We set w = v — v, y = dtw, and a = /j, — ]2. The 
function y is a solution of: 



d t y = DAy + fiy + ad t v in ^ to , 

-iti 

2/(0, ar) = crui(x) in Q, 



y(t,x) = on (7) 



The function rj(x, t) attains its minimum value with respect to the time at 
t = T' = ta 1 2 tl ■ We set tp = e~ sn y. Using the operator M 2 , introduced in 
Lemma [6TTT we introduce, following [23] (see also [24] and [25]). 

1 = J [ M 2 ip V- 
Let Ao be fixed as in Theorem 16.1 



t Jn 
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Lemma 6.2 Let A > Xq. There exists a constant C such that 

f"J 

Proof: From the Holder inequality, we have 



\1\<C 



S 3/2 A 2 / / e -a«, 3 3 + s -3/2 A -2 



2sr, 2 



<x 2 (W 



e" 2s y 



\1\ < s- 3/2 A- 2 f / T / (M 2 ^) 2 ) ( s 3 X 4 C 
\Jo Jn J y Jo 

Thus using Young's inequality, we obtain 



m < 7^ 3/2 a- 2 ( iim 2 ^ii 2 2(q , 1)+s 3 a 4 1 1 e - 2 *v 



1/2 



y 



(8) 



(9) 



Applying inequality a) of Lemma 16.11 to q := y, we obtain that there exists 
C > 0, such that 



„3\4 f*i 
'to 



< c 



*0 



(10) 



Furthermore, since fj, is bounded, and since ip is bounded from below by a 
positive constant, independent of A, we get that 



2e~ 2s W < « 3 A 4 / e- 2sl W 



(11) 



(12) 



'Qt "•tQ 

for A large enough. Combining (|10[) and (HJ), we obtain: 

IIM^II^+^A 4 /^ e - 2s VV 
< C [ S 3 A 4 j£ £ e~ 2 >y + / f Qll e-^, ^ {dt yf 

The conclusion of Lemma 16.21 follows from @ and (fT2"]) . □ 
Lemma 6.3 Let A > Ao. There exists a constant C such that 

+ A - 1 A- 2 e~ 2s ^ T '^a 2 {d t v) 2 
+ Jo e- 28 '"^'^) (DA w(T, x) + ^w{T' ', rr)) 2 ^ 

Proof: Using integration by parts over f2 and the boundary condition ip = 
on S*J , we get: 

—sXD L T , V ■ (^V/3)^ 2 + 2sA 2 L> J QT , <p\ V/3| 2 V> 2 . 
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We then obtain 



- [ ip (T', -f = X - s\ 2 D [ ip\VP\ 2 ip 2 dxdt + sXD [ tpAflip 2 dxdt 
2 Jn ■>< J QT ' 

since ip(to) = and Vip = XipWf3. As a consequence, for A > 1, and since A/3 is 
bounded in ft, we finally get 



-2srj{T' ,x) frpl n2 



y(T',x) 2 dx < 2\J\ +Cs\ 2 



T' 



t Jn 



(14) 



for some constant C. Using tp < ^ ^ <p 3 and Lemma [6TTT we get that: 



sX 2 



" 2s W < C 



A 



t Jul 



e- 2s \ 3 y 2 + s- 2 \- 2 / e- 2sr > 2[fx 2 y 2 + a 2 (d t vf 

J Ql 

(15) 

Arguing as in the proof of Lemma 16721 for equation (lll| . and since, for all x G ft, 
the function t t— > rj(t,x) attains its minimum over (to,ti) at t = T", we finally 
obtain that, for A large enough, the last term in (| 14[) is bounded from above by: 



C 



A 



2s Vy 2 + s- 2 \- 2 



- 2s ^ T '^a 2 {d t vf 



t J u 

If we now observe that 

y (T', x) = DAw (T 1 , x) + pw (T', x) + ad t v (T', x) , 

we get: 

a 2 v (T\ xf < 2y (T', xf + 2 [DAw (T', x) + pw (T 1 , x)f , 
and, since s > 1, the estimate of lemma 15751 follows. □ 
Lemma 6.4 We have < v, v < ule Mtl , and \d t v\ < (D + M)ule Mtl in 



Proof: From the parabolic maximum principle, we know that v,v > in Qt - 
Let h be the solution of the ordinary differential equation 



h' = hM on R+, 
h(0) = m. 



(16) 



The function h is increasing and H(t, •) := h(t) is a supersolution of the equa- 
tions satisfied by v and v. As a consequence of the parabolic maximum principle, 
we have, 

< v,v < h(ti) = u7e Mt \ in Q\\. (17) 



Let us set p :— dtv. The function p satisfies: 

d t p = DAp + pp in Q, 

p(t, x) = on E, 

p(0, a;) = DAui(x) + pui in ft. 



(18) 
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Since u, 6 V, p(0,x) e L°°(fi). Moreover, Hi(t,x) := -(D + M)ule Mt and 
H 2 {t,x) := (D + M)ule Mt are respectively sub- and supersolutions of (fT8|) . The 
parabolic maximum principle leads to the inequalities, 

-(D + M)u~e Mtl <p<(D + M)u- e Mtl in Q**. (19) 

□ 

From Lemmas 16.31 and 16.41 and since a vanishes outside fix, it follows that 



ifix e- 2s "( T '^)a 2 (%,T') 2 - Cs-^-^D + M) 2 ufe 2Mtl ) dx < 

Cs 3 / 2 X 2 J* 1 J u e" 2s V 3 y 2 + C J n e - 2s ^ T '^ [(Aw(T', x)) 2 + w(T', x) 2 ] dx. 

(20) 

Lemma 6.5 The ratio Uj/ul being fixed, there exists r > 0, independent of ft, 
Ui, Ui and u[, such that v (X", •) > ulr in fix- 

Proof: Let £o < 1 be a smooth function in fi, such that £o = 1 in i3j and 
£o = in fi\S E . Let £ be the solution of 

d t t, = DA£ - M£ in Q, 

x) = on E, (21) 
£(0,ar) = {yd/ui)io(x) in fi. 

Let us set r := inf £(T ,x). From the strong parabolic maximum principle, 

e > o in and therefore, we get that r > 0, since fix is a closed subset of 
fi. Moreover, the parabolic maximum principle also yields ul£ < v in In 
particular, we get v (T\ x) > ulr > in fix. □ 

p2C ( D + M)e Mtl 

From Lemma 16.51 it follows that, for A > \j , the term 

V s r 

v(x,T') 2 - Cs^X^iD + M) 2 m 2 e 2Mtl in flU]) satisfies: 

v(x,T') 2 - Cs- 1 \- 2 {D + M) 2 ufe 2Mtl > uf— > in fi x . 

We deduce that, for A large enough, 

In, e- 2sr,{T '' x) (J 2 dx 

< ^ ( S 3/2 A 2 Jtt ^ e -2s V(fi 3 y 2 + e -2« J (T',x)^ Au; ( T / i ^2 + ^jv ^2)^ 

(22) 

Using the fact that e 2sr] <p 3 remains bounded in QA we finally obtain 

Iklll^no < |s + l n ((Aw(T\x)) 2 + w (T',x) 2 ) dx^j . (23) 

Recalling that w = v — v, y = dt (v — v) , (|23[) implies the result of Theorem 



17 



6.2 Appendix B: Proof of Theorem 12.41 

Let it 7 be the solution of (Pn 7 ), and let v be the solution of (P^o). Let us set 



z 7 := v — Uy. The function z 7 is a solution of 



dtz 1 — DAz 7 — Zyfi(x) + 7^ 7 in Q, 
z 7 (t,i)=0onE, (24) 
z 7 (0, x) = in fl. 

It follows from the parabolic maximum principle that z^(t,x) > in Q\\. Thus 
u-y < v in Qll- Using the result of Lemma 16.41 we thus obtain: 

v»y(t,x) < ule Mtl in Q*J. (25) 

Let fc be the solution of 

k' = kM + 7W 2 e 2Mtl on R + , 
k(0) = 0. [2b > 

Since K(t,x) := kit) is a supersolution of (HM)) , we obtain that z 7 (f, x) < /c(£) 
in Qq 1 • Thus, since k is increasing, 

^^7-2 2 Aft i 

*r(t) < fc(ti) = 7 * M (e Mh 1) in Q*J. (27) 

Standard parabolic estimates (see e.g. [TB]) then imply, using Q17p. (|25l) . 
(j2~?| and the hypothesis ||M,|| C2 (jj) < Wj that: 



ll 9tV lli 2 (Q*J) = (llM«ll i 2 (Q *i ) + ||wi|| H i (n ))=0(W), (28) 
ll^t«rllia(0*j) = C , (ll^7A 1 + 7^ll L 2 (Q n ) ) = C^ 2 ), (29) 

H^lli^Q'i) = (II"t(M-7W7)II^(q«1) + lk||Hi(n)) 

= 0{ut), (30) 
sup \\v(t,.)\\ H 2 (n) = O(\\fiv\\ H i {0M . L 2 (n)) + \\ui\\ H 2^ n) ) = 0(ui) (31) 

*0<*<*1 

sup \\z 7 (t,.)\\ H 2 (n) = 0(\\z 7 n + ju 2 \\ H i {Q ^ tl . L 2 (n)) ), 
t <t<ti 

= 0(uf). (32) 

Using dHJ) and 122D, we get: 

G^,(i) = 0(u- 4 ). (33) 

Moreover, since uq = v, we have 

G>(7jA) = G>(G,/2) + G m (7,m) + 2(d t Zy,d t v - d t v) L2 ^ Q H ) 

+2(Az 1 (T',-),Av(T',.)-Av(T',-)) L2{n) "" (34) 
+2(^(T',-)^(T',-)-^(T" ) -))^(o)- 

Thus, using H§ and (J3TJ) , which are true for both u and u, and (gSJ), (021), ©, 
together with Cauchy-Schwarz inequality, we get: 

|G /1 ( 7 ,/i)-G' M (0,/i)| =0(?H 3 ). (35) 

□ 
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6.3 Appendix C: Definitions of the state spaces E and of 
their neighbourhood systems. 

In examples 1 and 2 of section 4, E is denned by 



where Xfc are the characteristic functions of the intervals (j§q, and a k are 

real numbers taken in finite subsets of [to, M]. 

In example 1, G { — 1,2}, and two distinct elements v\, v 2 of E, with z^i = 
YlLio a i,kXk(x) and ^ 2 = Sfc=io a 2,kXk(x) in S7i, are defined as neighbours if 
and only if there exists a unique integer ko in [10,89] such that ai : k 4 a 2,k - 
Note that, in this case, the number of elements in E is 2 80 . 

In example 2, a/. E {(M — m)j/20 + to, for j = . . . 20.}, and the neigh- 
bourhood system is defined as follows: two distinct elements V\, v 2 of E, with 

= Sfe=io a i,fcXfc(a;) and zy 2 = Y?k=w a 2,kXk(x) in Oi, are neighbours if and 
only if (i) there exists a unique integer ko in [10, 89] such that a\_k a 4 ot 2 , fc i (ii) 
additionally, \ai t k — ot 2 ,k \ = (M — m)/20. Note that, in such a situation, the 
number of elements in E is 21 80 . 

In examples 3 and 4 of section 4, E is defined by 

f " 1 

£7 := <p£ A-l, p(a:) = ^ aijXiji x ) m and p(a;) = to in fi\fli > , 

I «= 2 J 

where Xi,j are the characteristic functions of the square cells (i/20, i/20+1/20) x 
(j/20, j'/20+ 1/20), and ctij are real numbers taken in finite subsets of [to, M]. 

In example 3, Qjj E { — 1,2}. In this case, the number of elements in E 
is 2 256 . Two distinct elements v\, v 2 of E, with v\ = ^2ij=2 a i,jXi,j{ x ) and 

v 2 = Si,j=2 a i,jXi,j( x ) f° r £ E fii, are defined as neighbours if and only if there 
exists a unique couple (io, jo) of integers comprised between 2 and 17 such that 
a 1 -4 a 2 

to, JO I io,Jo 

In example 4 a,i^ E {(M — m)j/20 + to, for j = . . . 20.}. The number of 
elements in E is 21 256 . In this case, two distinct elements v\, v 2 of E, with 

v i = J2ij=2 a i,jXi,j(x) and v 2 = ^2]j =2 a i,jXi,j( x ) for x e fi i> are defined 
as neighbours if and only if (i) it exists a unique couple (io,jo) of integers 
comprised between 2 and 17 such that a£ - 7^ a ? ,j i and (ii) additionally 



K,,o-«Lol = ( M -^)/ 20 - 
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